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Abstract 

Certain families of combinatorial objects admit recursive descriptions in terms of gen- 
erating trees: each node of the tree corresponds to an object, and the branch leading to 
the node encodes the choices made in the construction of the object. Generating trees lead 
to a fast computation of enumeration sequences (sometimes, to explicit formulae as well) 
and provide efficient random generation algorithms. We investigate the links between 
the structural properties of the rewriting rules defining such trees and the rationality, 
algebraicity, or transcendence of the corresponding generating function. 

1 Introduction 

Only the simplest combinatorial structures — like binary strings, permutations, or pure invo- 
lutions (i.e., involutions with no hxed point) — admit product decompositions. In that case, 
the set Cl n of objects of size n is isomorphic to a product set: J7 n = [1, ei] x [1, e<2\ x • • • x [1, e n ]. 
Two properties result from such a strong decomposability property: (i) enumeration is easy, 
since the cardinality of Q n is e\C2 ■ ■ ■ e n ; (ii) random generation is efficient since it reduces to a 
sequence of random independent draws from intervals. A simple infinite tree, called a uniform 
generating tree is determined by the ef. the root has degree ei, each of its e± descendants has 
degree e2, and so on. This tree describes the sequence of all possible choices and the objects 
of size n are then in natural correspondence with the branches of length n, or equivalently 
with the nodes of generation n in the tree. The generating tree is thus fully described by its 
root degree (ei) and by rewriting rules, here of the special form, 

(ei) ^ (e i+ i) (e i+ i) ■ ■ ■ (e i+1 ) = {e i+ if\ 

where the power notation is used to express repetitions. For instance binary strings, permu- 
tations, and pure involutions are determined by 

S : [(2), (2) - (2) (2)] 

V : [(1), {(fcW(fc + l)*}*>i] 

J : [(1), {(2k-l)^(2k + l) 2k - 1 } k > 1 ]. 



"Corresponding author. 



1 



A powerful generalization of this idea consists in considering unconstrained generating 
trees where any set of rules 

s = K s o), {(fc) ~» (ei,fc) (e 2 ,fc) • • • (efe,fc)}] (1) 

is allowed. Here, the axiom (sq) specifies the degree of the root, while the productions e^k 
list the degrees of the k descendants of a node labeled k. Following Barcucci, Del Lungo, 
Pergola and Pinzani, we call S an ECO-system (ECO stands for "Enumerating Combinatorial 
Objects"). Obviously, much more leeway is available and there is hope to describe a much 
wider class of structures than those corresponding to product forms and uniform generating 
trees. 

The idea of generating trees has surfaced occasionally in the literature. West introduced it 
in the context of enumeration of permutations with forbidden subsequences |27| I28j: this idea 
has been further exploited in closely related problems [51 151 112| I13| . A major contribution in 
this area is due to Barcucci, Del Lungo, Pergola, and Pinzani [3J |H] who showed that a fairly 
large number of classical combinatorial structures can be described by generating trees. 

A form equivalent to generating trees is well worth noting at this stage. Consider the walks 
on the integer half-line that start at point (sq) and such that the only allowable transitions are 
those specified by £ (the steps corresponding to transitions with multiplicities being labeled). 
Then, the walks of length n are in bijective correspondence with the nodes of generation n in 
the tree. These walks are constrained by the consistency requirement of trees, namely, that 
the number of outgoing edges from point k on the half-line has to be exactly k. 

Example 1. 123-avoiding permutations 

The method of "local expansion" sometimes gives good results in the enumeration of per- 
mutations avoiding specified patterns. Consider for example the set S n (123) of permuta- 
tions of length n that avoid the pattern 123: there exist no integers i < j < k such that 
a(i) < a(j) < a(k). For instance, a = 4213 belongs to 64(123) but a = 1324 does not, as 
a(l) < a(3) < a(4). 

Observe that if r £ S n _|_i(123), then the permutation a obtained by erasing the entry 
n + 1 from r belongs to S n (123). Conversely, for every a £ (5 n (123), insert the value n + 1 in 
each place that gives an element of © n +i(123) (this is the local expansion). For example, the 
permutation a = 213 gives 4213, 2413 and 2143, by insertion of 4 in first, second and third 
place respectively. The permutation 2134, resulting from the insertion of 4 in the last place, 
does not belong to 64(123). This process can be described by a tree whose nodes are the 
permutations avoiding 123: the root is 1, and the children of any node a are the permutations 
derived as above. Figure ^a) presents the first four levels of this tree. 

Let us now label the nodes by their number of children: we obtain the tree of Figure ^b). 
It can be proved that the k children of any node labeled k are labeled respectively k + 
1,2, 3,... , k (see [2Zj). Thus the tree we have constructed is the generating tree obtained 
from the following rewriting rules: 

[(2), {(fc)^(2)(3)...(fc-l)(fc)(fc + l)} fc > 2 ]. 

The interpretation of this system in terms of paths implies that 123-avoiding permuta- 
tions are equinumerous with "walks with returns" on the half-line, themselves isomorphic 
to Lukasiewicz codes of plane trees (see, e.g., [2SJ p. 31-35]). We thus recover a classic 
result 123-avoiding permutations are counted by Catalan numbers; more precisely, 

|6 n (123)| = ( 2 ;)/(n + l). □ 
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4321 3241 4231 4213 2143 4312 3142 4132 
3421 3214 2431 2413 3412 1432 
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5234 32 423 42 3 32 




(a) 



(b) 



Figure 1: The generating tree of 123-avoiding permutations, (a) Nodes labeled by the per- 
mutations, (b) Nodes labeled by the numbers of children. 



We shall see below that (certain) generating trees correspond to enumeration sequences 
of relatively low computational complexity and provide fast random generation algorithms. 
Hence, there is an obvious interest in delineating as precisely as possible which combinatorial 
classes admit a generating tree specification. Generating functions condense structural infor- 
mation in a simple analytic entity. We can thus wonder what kind of generating function can 
be obtained through generating trees. More precisely, we study in this paper the connections 
between the structural properties of the rewriting rules and the algebraic properties of the 
corresponding generating function. 

We shall prove several conjectures that were presented to us by Pinzani and his coauthors 
in March 1998. Our main results can be roughly described as follows. 

- Rational systems. Systems satisfying strong regularity conditions lead to rational gen- 
erating functions (Section EJ). This covers systems that have a finite number of allowed 
degrees, as well as systems like (2. a), (2.6), (2.c) and (2.d) below where the labels are 
constant except for a fixed number of labels that depend linearly and uniformly on k. 

- Algebraic systems. Systems of a factorial form, i.e., where a finite modification of the 
set {1, ...,k} is reachable from k, lead to algebraic generating functions (Section |3J). 
This includes in particular cases (2./) and (2. g). 

- Transcendental systems. One possible reason for a system to give a transcendental series 
is the fact that its coefficients grow too fast, so that its radius of convergence is zero. 
This is the case for System (2.h) below. Transcendental generating functions are also 
associated with systems that are too "irregular". An example is System (2.e). We shall 
also discuss the holonomy of transcendental systems (Section 0J). 
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Example 2 . A zoo of rewriting systems 

Here is a list of examples recurring throughout this paper. 



[(3), {(k) ~> (3) fc - 3 (fc + l)(k + 2)(k + 9)}] (2.a) 

[(3), {(k) -(3)^(3^ + 6)}] (2.6) 

[(2), {(fc) ~> (2) fc ~ 2 (2 + (fc mod 2)){k + 1)}] (2.c) 

[(2), {(k) (2) fc - 2 (3 - (fc mod 2))(* + 1)}] (2.d) 

[(2), {(*:) - (2) fc - 2 (3 - [3p:fc = 2P])(k + 1)}] (2.c) 

[(2),{(fc)^(2)(3)...(fc-l)(fc)(fc + l)}] (2./) 

[(l),{(fc)^(l)(2)...(fc-l)(fc + l)}] (2. 5 ) 

[(2),{(A;)-(2)(3)(fe + 2) fc - 2 }] (2.h) 

(In (2.e), we make use of Iverson's brackets: [P] equals 1 if P is true, otherwise.) □ 



Notations. From now on, we adopt functional notations for rewriting rules: systems will 
be of the form 

[(so), {(k)^( ei (k))(e 2 (k))...(e k (k))}] 

where sq is a constant and each e; is a function of k. Moreover, we assume that all the values 
appearing in the generating tree are positive: each node has at least one descendant. 

In the generating tree, let f n be the number of nodes at level n and s n the sum of the 
labels of these nodes. By convention, the root is at level 0, so that /o = 1. In terms of walks, 
f n is the number of walks of length n. The generating function associated with the system is 

F(Z) = J2 fnZ n . 
n>0 

Remark that s n = f n +i, and that the sequence (f n ) n is nondecreasing. 

Now let f n ,k be the number of nodes at level n having label k (or the number of walks of 
length n ending at position k). The following generating functions will be also of interest: 

F{z,u) = J2 fn,kZ n u k and F k (z) = ^ /„/. 

n,fc>0 n>0 

We have F(z) = F(z, 1) = J2k>i Fk(z). Furthermore, the F^'s satisfy the relation 

F k (z) = [k = s ]+zJ2^,kF j (z), (2) 

where Wj jk = \{i < j : ej(j) = k}\ denotes the number of one-step transitions from j to k. 
This is equivalent to the following recurrence for the numbers 

fo,k = [k = s ] and f n +\,k = n j,kfn,j, (3) 
that results from tracing all the paths that lead to k in n + 1 steps. 
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Counting and random generation. The recurrence (j3J) gives rise to an algorithm that 
computes the successive rows of the matrix (f n ,k) by "forward propagation": to compute 
the (n + l)th row, propagate the contribution f n j to / n +i, ei (j) f° r an pairs (i,j) such that 

1 < j. Assume the system is linearly bounded: this means that the labels of the nodes that 
can be reached in m steps are bounded by a linear function of m. (All the systems given 
in Example EJ except for (2.6), are linearly bounded; more generally, systems where forward 
jumps are bounded by a constant are linearly bounded.) Clearly, the forward propagation 
algorithm provides a counting algorithm of arithmetic complexity that is at most cubic. 

For a linearly bounded system, uniform random generation can also be achieved in poly- 
nomial time, as shown in .2 ]. We present here the general principle. 

Let g n ^k be the number of walks of length n that start from label k. These numbers are 
determined by the recurrence g n ^ = J2i 9n-l,ei(k)i that traces all the possible continuations 
of a path given its initial step. Obviously, f n = g n .s , with so the axiom of the system. As 
above, the g n ^ can be determined in time 0(n 3 ) and 0(n 2 ) storage. Random generation 
is then achieved as follows: In order to generate a walk of length n starting from state k, 
pick up a transition i with probability g n -i,ei(k)/9n,k, and generate recursively a walk of 
length n — 1 starting from state ej(fc). The cost of a single random generation is then 0(n 2 ) if 
a sequential search is used over the O(n) possibilities of each of the n random drawings; the 
time complexity goes down to 0(n log n) if binary search is used, but at the expense of an 
increase in storage complexity of 0(n 3 ) (arising from 0(n 2 ) arrays of size 0(n) that binary 
search requires). 

2 Rational systems 

We give in this section three main criteria (and a variation on one of them) implying that the 
generating function of a given ECO-system is rational. 

Our first and simplest criterion applies to systems in which the functions ej are uniformly 
bounded. 

Proposition 1 If finitely many labels appear in the tree, then F{z) is rational. 

Proof. Only a finite number of Fk's are nonzero, and they are related by linear equations 
like Equation © above. a 



Example 3. The Fibonacci numbers 

The system [(1), {(k) ~> {k) k - l ({k mod 2) + l)}] can be also written as [(1), {(1) ~c (2), (2) ~» 
(1)(2)}]. Hence the only labels that occur in the tree are 1 and 2. Eq. © gives F\(z) = 
1 + zF 2 (z) and F 2 (z) = z{Fi{z) + F 2 {z)). Finally, 

F(z) = - 5 = V f n z n = 1 + z + 2z 2 + 3z 3 + 5z 4 + • • • , 

the well-known Fibonacci generating function. □ 

None of the systems of Example El satisfy the assumptions of Proposition ^ However, the 
following criterion can be applied to systems (2. a) and (2.6). 
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Proposition 2 Let o~(k) = e\(k) + e2(k) + • • • + et(k). If a is an affine function of k, say 
o~(k) = ak + P, then the series F(z) is rational. More precisely: 

1 + (sp - a)z 
[ ' ~ l-az-0z*' 

Proof. Let n > and let k\,k%, . . .k* denote the labels of the f n nodes at level n. Then 

fn+2 = Sn+i = (ah + (3) + (ak 2 +/3) H h (akf n + 0) 

= as n + (3f n = a/n+i + 0fn- 

We know that /o = 1 and f\ = sq. The result follows. a 



Example 4. Bisection of Fibonacci sequence 

The system [(2), {(fc) ~» (2) fe - 1 (A: + 1)}] gives F(z) = p^fp = 1 + 2z + 5z 2 + • • •, the 
generating function for Fibonacci numbers of even index. (Changing the axiom to (sq) = (3) 
leads to the other half of the Fibonacci sequence.) Some other systems, like 

[(2),{(k)^(l) k - l (2k)}}, 

[(2), {(Jfe) -s* (2) fc ^ 2 (3 - (Jfe mod 2))(fc + (k mod 2))}], 
[(2), {{k) ~» (2) fc ~ 2 (3 - [jfe is prime]) (A; + [jfc is prime])}], 

lead to the same function F(z) since cr(fc) = 3fc — 1 and sq = 2. However, the generating trees 
are different, as are the bivariate functions F(z,u). □ 



Example 5. Prime numbers and rational generating functions 

Amazingly, it is possible to construct a generating tree whose set of labels is the set of prime 
numbers but that has a rational generating function F(z). This is a bit unexpected, as 
prime numbers are usually thought "too irregular" to be associated with rational generating 
functions. For n > 1, let p n denote the nth prime; hence (pi,P2,P3, ■ ■ ■) = (2,3,5,...). 
Assume for the moment that the Goldbach conjecture is true: every even number larger than 
3 is the sum of two primes. Remember that, according to Bertrand's postulate, p n +i < 2p n 
for all n (see, e.g., [22 P- 140]). 

For n > 1, the number 2p n — p n +\ + 3 is an even number larger than 3. Let q n and r n 
be two primes such that 2p n — p n +i + 3 = q n + r n . In particular, q% = r\ = 2. Consider the 
system 

[(2),{(Pn)-(Pn+l)(9n)(r„)(2r- 3 }]. 

It satisfies the criterion of Proposition [21 with o~{k) = Ak — 3. Hence, the generating function 
of the associated generating tree is 



l-2z 1 
{Z) ~ 1 - 4z + 3z 2 ~ 2 



1 1 

+ 



1-z 1-3* 



Consequently, the number of nodes at level n is simply f n = (1 + 3 n )/2. This can be checked 
on the first few levels of the tree drawn in Figure El 

Now, one can object that the Goldbach conjecture is not proved; however, it is known 
that every even number is the sum of at most six primes j^lj, and a similar example can be 
constructed using this result. 

□ 
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Figure 2: A generating tree with prime labels and rational generating function. 

Proposition |2] can be adapted to apply to systems that "almost" satisfy the criterion of 
Proposition |21 like System (2.c) or (2.d). Let us consider a system of the form 

(so), (k) ~> e\\k), . . . , e^\k) if k is even, 
(k)^> e l ? ] (k),...,e% ] (k) if k is odd. 

Assume, moreover, that: 

(i) the corresponding functions <ro and o~\ are afhne and have the same leading coefficient 
a, say a^(k) = ak + [3 and ai(k) = ak + j3\; 

(ii) exactly m odd labels occur in the right-hand side of each rule, for some m > 0. 

Proposition 3 If a system satisfies properties (i) and (ii) above, then 

p.. _ l + (s - a)z + (s\ - asp - (3 )z 2 
1 — az — Pqz 2 — m((3\ — 0o)z 3 

Of course, if [3$ = f3\ , we recover the generating function of Proposition |21 

Proof. The proof is similar to that of Proposition |^1 The only new ingredient is the fact 
that, for n > 1, the number of nodes of odd label at level n is mf n -\. ■ 

System (2.c) satisfies properties (i) and (ii) above with a = 3, (5q = —1, /?i = 0, m = 1, 
so = 2 and s\ = 5. Consequently, its generating function is F(z) = 1 _ 3 ^^_ z n ■ System (2.d), 
although very close to (2.c), does not satisfy property (ii) above, so that Proposition |31 does 
not apply. However, another minor variation on the argument of Proposition based on the 
fact that the number o n of odd labels at level n satisfies o n = 2(/ n _i — o n _i), proves the 
rationality of F(z). 

Alternatively, rationality follows from the last criterion of this section, which is of a 
different nature. We consider systems [(sq), {(k) (e±(k))(e2(k)) . . . (ek(k))}] that can be 
written as 

[(so), {(A;) ~» (ci(fe))(p2(A;)) . . . (ck- m (k))(k + ai )(k + a 2 )...(k + a m )}\ (4) 

where 1 < a\ < a 2 < • • • < a m and the functions Cj are uniformly bounded. Let C = 
m&x ijk {s ,Ci(k)}. 
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Proposition 4 Consider the system Ijl]). and let irj tk = \{i < j : ej(j) = fc}|. If all th 



e series 



J2 *j,k t J 

for k < C are rational, then so is the series F(z). 

Proof. We form an infinite system of equations defining the series F k {z) by writing Eq. J2J) 
for all k > 1. In particular, for k > C, we obtain 

m 

F k (z) = zY,F k - at {z), 

i=\ 

with Fj(z) = if j < 0. This part of the system is easy to solve in terms of F±, . . . ,Fq- 
Indeed, for fcgZ: 

c 

F k (z) = ^2P i , k (z)F i (z) (5) 
i=i 

where the Pi k are polynomials in z defined by the following recurrence: for all i < C, 





[fe = i] 



if fe < 0, 

if < k < C, 



ni 



zJ2 Pi,k-a e (z) Hk>C. 



(6) 



Using ©, we find 



c 



F(z) = Y,F k (z)=Y: 



k>l 



4 = 1 



k>l 



According to ©, for all i < C, the series X)fc>i Pi,k( z )t k is a rational function of z and t, of 
denominator 1 — z ^£ • At t = 1, it is rational in z. Hence, to prove the rationality of F(z), 
it suffices to prove the rationality of the Fi(z), for i < C. 

Let us go back to the C first equations of our system; using @, we find, for k < C: 



c 



F k (z) = [k = s ] + z^2 



i=l 



Fi(z)^P iyj (z)TT^ k 



Again, J2j>i Pi,j( z ) 7T j,kt J 1S a rational function of z and t (the Hadamard product of two 
rational series is rational). Thus the series F k {z), for k < C, satisfy a linear system with 
rational coefficients: they are rational themselves, as well as F(z). m 

Examples (2. a), (2.c), (2.d) and (2.e) have the form (0J). The above proposition implies 
that the first three have a rational generating function. System (2.e) will be discussed in 
Sectional and proved to have a transcendental generating function. 
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3 Factorial walks and algebraic systems 



In this section, we consider systems that are of a factorial form. By this, we mean informally 
that the set of successors of (k) is a finite modification of the integer interval {1,2, ... ,k}. 
As was detailed in the introduction, ECO-systems can be rephrased in terms of walks over 
the integer half-line. We thus consider the problem of enumerating walks over the integer 
half-line such that the set of allowed moves from point k is a finite modification of the 
integer interval [0,k]. We shall mostly study modifications around the point k (although 
some examples where the interval is modified around as well are given at the end of the 
section). Precisely, a factorial walk is defined by a finite (multi)set A C Z and a finite set 
B C N + , where N + = {1, 2, 3, . . .}, specifying respectively the allowed supplementary jumps 
(possibly labeled) and the forbidden backward jumps. In other words, the possible moves from 
k are given by the rule: 

(k) ^ [0, k - 1] \ (k - B) U (k + A). (7) 

Observe that these walk models are not necessarily ECO-systems, first because we allow labels 
to be zero - but a simple translation can take us back to a model with positive labels - and 
second because we do not require (A;) to have exactly k successors. 

We say that an ECO-system is factorial if a shift of indices transforms it into a factorial 
walk. Hence the rules of a factorial ECO-system are of the form 

(fc + r) ~> [r, k + r - 1] \ (k + r - B) U (k + r + A), 

that is, 

(k)^[r,k-l]\(k- B) U (k + A) for k > r > 1. (8) 

The generating function F(z) for such an ECO-system, taken with axiom (so), equals the 
generating function for the walk model (|7j). taken with axiom (sq — r). However, remember 
that the rewriting rules defining a generating tree have to obey the additional condition that 
a node labeled k has exactly k successors. Taking k = r in l|8"|). this implies that r = |^4|. 
Taking k > r + max 5, this implies that r + \B\ = \ A\, so that finally 5 = 0. Hence, strictly 
speaking, either one has a "fake" factorial ECO-system (that is some of its initial rules are 
not of the factorial type), either one has a "real" factorial ECO-system and then it is given 
by rules of the form 

(k) [r, k - 1] U (k + A) for k > r > 1, 

where A is a multiset of integers of cardinality r. For instance, Systems (2./) and (2. g) are 
factorial. We shall prove that all factorial walks have an algebraic generating function. The 
result naturally applies to factorial ECO-systems. 

We consider again the generating function F(z,u) = J2n k>o fn,kZ n u k , where f n ^ is the 
number of walks of length n ending at point k. We also denote by F^{z) the coefficient of 
u k in this series, and by f n (u) the coefficient of z n . The first ingredient of the proof is a 
linear operator M, acting on formal power series in u, that encodes the possible moves. More 
precisely, for all n > 0, we will have: 

M[f n ](u) = f n+1 (u). 
The operator M is constructed step by step as follows. 
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The set of moves from k to all the positions 0, 1, . . . , k — 1 is described by the operator 
Lq that maps u k to u° + u 1 + • • • + u k ~ 1 = (1 — u k ) /(l — u). As Lq is a linear operator, 
we have, for any series g{u): 

r , ,/ i g(l) -g(-g) 
L o[g}(u) = — — - — . 

The fact that transitions near k are modified, with those of type k + a (with a £ A) 
allowed and those of type k — (3 (with G B) forbidden, is expressed by a Laurent 
polynomial 



P(u) = PkU k = A(u) - B(u) with A(u) = ^u a and B(u) = ^ u 



/3 . 

k=-b a£A /3eB 



The degree of P is a := max A, the largest forward jump and b := max(0, — B, — A) 
is largest forbidden backward jump or the largest supplementary backward jumps (we 
take b = if the set B is empty) . 

The operator 

L[g}(u) := L [g}(u) + P{u)g{u) 
describes the extension of a walk by one step. 
— Finally, the operator M is given by 

M[g}(u)=L[g}(u)-{u<°}L[g}(u), 

where {u <0 }h(u) is the sum of all the monomials in h(u) having a negative exponent. 
Hence M is nothing but L stripped of the negative exponent monomials, which corre- 
spond to walks ending on the nonpositive half-line. Observe that, for any series g(u), 
the only part of L[g}{u) that is likely to contain monomials with negative exponents is 
P(u)g(u). Consequently, 

M[g](u) = L[g](u) - {u<°}[P(u)g(u)] 
and if g(u) = J2k9kU k , then 

{u< }[P(u)g(u)} =Y / T / 9kP- i u k - 1 = "£g k r k (u). (9) 

i=l k=0 k=0 

Assume for simplicity that the initial point of the walk is 0; other cases follow the same 
argument. The linear relation / n +i(n) = M[f n ](u), together with fo(u) = 1, yields 

F(z,u) = 1 + zM[F](z,u) 

= i + z 1) ~ F ^ u ) + P( u )F(z, u) + {u <0 }[P(u)F(z, u)\ 

V 1 — u 



Thanks to ©, we can write 



6-1 

{u< }[P(u)F(z,u)] = J2n(u)F k (z 

k=0 
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where rfc(ii) is a Laurent polynomials (defined by Equation |!5J) whose exponents belong to 
[k — b, — 1]. Thus, F(z,u) satisfies the following functional equation: 

F(z, „) (l + - zP(u)) = 1 + ^E^L + z £ r k {u)F k {z). (10) 

\ II 1 II J 



Let us take an example. The moves 

(k) ~» (0)(1) ■■■ (k — 5)(fc - 3)(fc - l)(fc)(fc + 7) (A; + 9), 
lead to -A(u) = u + u 7 + u 9 and -B(u) = u™ 4 + u~ 2 . Moreover, 

{u <0 }[B(u)F(z, u)] = (u~ 2 + u~ 4 )F (z) + (n" 1 + u -3 )^^) + u~ 2 F 2 {z) + u" 1 ^), 
so that the functional equation defining F(z, u) is 

F(z, u) (l + YZT^ - z (! + u 7 + u 9 - u~ 4 - U - 2 )^j = 

l + ^1) + z(n _ 2 + n -4 )i?()(z) + + u -3 )Fl(z) + ZU -2 F2 ( Z ) + zu -lF 3 (z). 



The second ingredient of the proof, sometimes called the kernel method, seems to belong 
to the "mathematical folklore" since the 1970's. It has been used in various combinatorial 
problems |1U1 1181 120j and in probabilities ^3] . See also jSJ for more recent and system- 
atic applications. This method consists in cancelling the left-hand side of the fundamental 
functional equation ()1U|) by coupling z and u, so that the coefficient of the (unknown) quantity 
F(z, u) is zero. This constraint defines u as one of the branches of an algebraic function of z. 
Each branch that can be substituted analytically into the functional equation yields a linear 
relation between the unknown series F{z, 1) and F^z), < k < b. If enough branches can be 
substituted analytically, we obtain a system of linear equations, whose solution gives F(z, 1) 
and the F k (z) as algebraic functions. From there, an expression for F(z,u) also results in the 
form of a bivariate algebraic function. 

Let us multiply Eq. (|10j) by u b (l — u) to obtain an equation with polynomial coefficients 
(remind that we take b = if the set B of forbidden backward steps is empty). The new 
equation reads K(z,u)F(z,u) = R(z,u), where K(z,u) is the kernel of the equation: 

K{z,u) = u b (l - u) (l + - zP(u)^j , 

= u^l-^+z^-zil-^Y.^ + zil-u)^^'' 3 - (11) 

This polynomial has degree a + b + 1 in u, and hence, admits a + b + 1 solutions, which are 
algebraic functions of z. The classical theory of algebraic functions and the Newton polygon 
construction enable us to expand the solutions near any point as Puiseux series (that is, series 
involving fractional exponents; see JJ)- The a + b + 1 solutions, expanded around 0, can be 
classified as follows: 

- the "unit" branch, denoted by uq, is a power series in z with constant term 1; 
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- 6 "small" branches, denoted by u\, . . . , «&, are power series in z l / b whose first nonzero 
term is (z 1 ^, with ( b + 1 = 0; 

- a "large" branches, denoted by v\, . . . , v a , are Laurent series in z 1//a whose first nonzero 
term is £z -1 / a , with ( a + 1 = 0. 

In particular, all the roots are distinct. (It is not difficult to check "by hand" the existence 
of these solutions: for instance, plugging z = t b and u = tw(t) in K(z,u) = confirms the 
existence of the 6 small branches.) Note that there are exactly 6 + 1 finite branches: the unit 
branch uq and the b small branches m, . . . , As F(z, u) is a series in z with polynomial 
coefficients in u, these 6+1 series Uj, having no negative exponents, can be substituted for 
u in F(z,u). More specifically, let us replace u by Ui in (|1()[) : the right-hand side of the 
equation vanishes, giving a linear equation relating the 6+1 unknown series F(z, 1) and 
Fk(z), < k < b. Hence the 6+1 finite branches give a set of 6 + 1 linear equations relating 
the 6 + 1 unknown series. One could solve directly this system, but the following argument is 
more elegant. 

The right-hand side of ()1U|) . once multiplied by u b (l — u), is 



R(z,u) = u b {l -«) 1 + 



6-1 \ 

F(z,l) + zJ2r k (u)F k (z) . 



By construction, it is a polynomial in u of degree 6 + 1 and leading coefficient —1. Hence, it 
admits 6 + 1 roots, which depend on z. Replacing u by the series uo, iti, . . . , Ub in Eq. (jlOf) 
shows that these series are exactly the 6 + 1 roots of R, so that 

6 

R(z,u) = - Y[(u - Ui). 

i=0 

Let p a := [u a ]P(u) be the multiplicity of the largest forward jump. Then the coefficient of 
u a+b+1 in K(z, u) is p a z, and we can write 

6 a 

K(z,u) = p a zY[(u - Ui)Y[(u - Vi). 

1=0 8=1 

Finally, as K(z,u)F(z,u) = R(z,u), we obtain 

F(z u) = -Ut=o(u-u t ) = 1 

V ' ' u b (l-u)+zu b -zu b (l-u)P(u) P a zlli=i(u-Vi)' 

We have thus proved the following result. 

Proposition 5 The generating function F(z, u) for factorial walks defined by (J7J) and starting 
from is algebraic; it is given by (|12j) . where uq,. . . , (resp. vi, . . . , v a ) are the finite (resp. 
infinite) solutions at z = of the equation K(z,u) = and the kernel K is defined by (|llj) . 
In particular, the generating function for all walks, irrespective of their endpoint, is 



1 b 



^i) = -tI1( 1 



Ui 



Z n 
1=0 
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and the generating function for excursions, i.e., walks ending at 0, is, for b < 0: 



b b 



Z " 

1=0 

(for 6 = 0, the relation becomes F(z,0) = i+ z ^ P(jZ T\i=o u i-) 

These results could be derived by a detour via multivariate linear recurrences, and the present 
treatment is closely related to (Hll^; however, our results were obtained independently in 
March 1998 Q]. 

The asymptotic behaviour of the number of n-step walks can be established via singularity 
analysis or saddle point methods. The series Ui have "in general" a square root singularity, 
yielding an asymptotic behaviour of the form A[i n n~ 3 / 2 . We plan to develop this study in a 
forthcoming paper. 

Example 6. Catalan numbers 

This is the simplest factorial walk, (k) (0)(1) . . . (k)(k + 1), which corresponds to the 
ECO-system (2.f). The operator M is given by 

M[f]{ u ) = WLJM + (i + u) f {u ). 

1 — u 

The kernel is K(z, u) = 1 — u + z — z(l — u)(l + u) = 1 — u + zu 2 , hence Uq(z) = 1- ^~ 42 , so 
that 

T?( -n l-np _ \ -2z- y/T=Az~ _ ^ (2n\ z n ~ l 
[Z > } ~ z ~ 2z 2 "^lnln + 1' 

n>l \ / 

the generating function of the Catalan numbers (sequence M1459 1 ). This result could be 
expected, given the obvious relation between these walks and Lukasiewicz codes. □ 



Example 7. Motzkin numbers 

This example, due to Pinzani and his co-authors, is derived from the previous one by forbid- 
ding "forward" jumps of length zero. The rule is then 

(k) ~> (0) • • • (A; — l)(k + 1). 

The operator M is 

The kernel is K(z, u) = 1 — u + z — zu(l — u) = 1 + z — u(l + z) + zu 2 , leading to 



F(Z,1) = 1 g VI ^ 2 ^ f! = 1 + z + 2z 2 +4z 3 + 9z 4 + 21z 5 + 0(z 6 ); 

the generating function for Motzkin numbers (sequence Ml 184). □ 



lr The numbers Mxxxx are identifiers of the sequences in The Encyclopedia of Integer Sequences |24| 
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Example 8. Schroder numbers 

This example is also due to the Florentine group. The rule is (k) ~> (0) . . . (k— l)(k)(k + l) 2 . 
From Proposition [5J we derive 



F ( Zj 1) = 1 3Z ^ fe _±i! = 1 + 3z + H Z 2 + 45z 3 + 197z 4 + ^5) 

The coefficients are the Schroder numbers (M2898: Schroder's second problem). We give 
in Table 1 at the end of the paper a generalization of Catalan and Schroder numbers, corre- 
sponding to the rule (k) ~> (0) . . . (k — l)(k)(k + l) m . This generalized rule has recently been 
shown to describe a set of permutations avoiding certain patterns |19j . □ 

The above examples were all quadratic. However, it is clear from our treatment that 
algebraic functions of arbitrary degree can be obtained: it suffices that the set of "exceptions" 
around k have a span greater than 1. Let us start with a family of ECO-systems where 
supplementary forward jumps of length larger than one are allowed. 

Example 9. Ternary trees, dissections of a polygon, and m-ary trees 
The ECO-system with axiom (sq) = (3) and rule 

(k) ~> (3) (4) • • • (k)(k + l)(k + 2) 

is equivalent to the walk 

(k) ~> (0)(1) • • • (k)(k + l)(k + 2). 
The kernel is K(z, u) = 1 — u + zu 3 , and the generating function 

/3n\ z n ~ l 

n>l ' 



n I 2n + 1 



counts ternary trees (M2926). 

More generally, the system with axiom (m) and rewriting rules 

(jfe) ~» (m) • • • (*;)(*! + l)(Jfe + 2) • • • (jfe + m - 1) 

yields the m-Catalan numbers, — l)n + 1), that count m-ary trees. The kernel 

is 1 — u + zu m and the generating function F(z, 1) satisfies F(z, 1) = (1 + zF(z, l)) rn . In 
particular, the 4-Catalan numbers ( 4 ^)/(3n + 1) appear in .24j (sequence M3587) and count 
dissections of a polygon. 

□ 

In the above examples, all backward jumps are allowed. In other words, each of these 
examples corresponds to an ECO-system. Let us now give an example where backward 
jumps of length 1 are forbidden. 

Example 10. 

Consider the following modification of the Motzkin rule: 

(k) ~> (())■■■ (jfe- 2)(lfe + l). 
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The kernel is now K(z,u) = u(l — u) + zu — z(l — u)(u 2 — 1), and, according to 1)12(1 . the 
series F(z) = F(z, 1) is given by F(z) = l/[z(vi — 1)], where v\ satisfies K(z,v±) = and is 
infinite at z = 0. Denoting G = zF(z), we find that the algebraic equation defining G is: 

1 + 2G + G 2 + G 3 

G = z TTg ' 

□ 

So far, we have only dealt with walks for which the set of allowed moves was obtained by 
modifying the interval [0, k] around k. One can also modify this interval around 0: we shall 
see - in examples - that the generating function remains algebraic. However, it is interesting 
to note that in these examples, the kernel method does not immediately provide enough 
equations between the "unknown functions" to solve the functional equation. 

Let us first explain how we modify the interval [0, k] around 0. The walks we wish to count 
are still specified by a multiset A of allowed supplementary jumps and a set B of forbidden 
backward jumps. But, in addition, we forbid backward jumps to end up in C, where C is a 
given finite subset of N. In other words, the possible moves from k are given by the rule 

(k) ~» [0, k - 1] \ (C U (k - B)) U (k + A). 

Again, we can write a functional equation defining F(z,u): 

F(z,u) = l + z ( nzA) -F(z,u) + p{u)F{z ^ u) + g r k {u)F k {z) - ]T u^G^z) ] , (13) 



where, as above, 



p(u) = u a - u ~ p and niu) = Y uk ~^ 

aeA j3>k, /3eB 

the new terms in the equations being 

G 7 (z) = F(z, 1) - £ F k (z) - £ F/3 +J (z). 

k=0 /3eB 

Observe that the first three terms are the same as in the case C = 0. The equation, once 
multiplied by u b (l — u), reads K(z,u)F(z,u) = R(z,u) where K(z,u) is given by (fTTj) and 



R(z, u) = u\l -u) [1 + + z Y r k (u)F k (; 




fc=0 

The kernel is not modified by the introduction of C. As above, it has degree a + 6+1 in u, and 
admits b + 1 finite roots uq, . . . , Uf, around z = 0. However, R(z, u) now involves b + 1 + \C\ 
unknown functions, namely F{z, 1), the F k (z), < k < b and the Gj(z), 7 S C. The degree 
of R in u is no longer 6 + 1 but b + c + 1, where c = maxC. The 6 + 1 roots of K that 
can be substituted for u in Eq. (fT3)) provide 6 + 1 linear equations between the 6 + |C| + 1 
unknown functions. Additional equations will be obtained by extracting the coefficient of vP 
from Eq. (|13|). for some values of j. In general, we have: 

F 3 (z) = \j = 0] + z Y Fj-a(z) + z[j ? C] I F(z, 1) - £ F k (z) - ^ F j+P {z) I . (14) 
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It is possible to construct a finite subset S C N such that the combination of the b + 1 
equations obtained via the kernel method and the equations ()14p written for j £ S determines 
all unknown functions as algebraic functions of z - more precisely, as rational functions of z 
and the roots Uq, . . . , uj, of the kernel. However, this is a long development, and these classes 
of walks play a marginal role in the context of ECO-systems. For these reasons, we shall 
merely give two examples. The details on the general procedure for constructing the set S 
can be found in [Jj. 

Example 11. 

This example is obtained by modifying the Motzkin rule of Example around the point 0. 
Take A = C = {1} and B = 0. The rewriting rule is 

(fc)^(0)(2)(3)---(fc-l)(fc + l). 

The functional equation reads 

(1 - u + z - zu(l - u))F{z, u) = l-u + zF(z,l) - zu(l - u)Gx(z), (15) 

with G\{z) = F(z, 1) — Fq(z) — Fi(z). The kernel has a unique finite root at z = 0: 

1 + z - Vl - 2z - 3z 2 



2z 



whereas the right-hand side of Eq. (|15|) contains two unknown functions. Writing Eq. ()14|) 
for j = and j = 1 yields 

F Q (z) = l + z(F(z,l)-F (z)) and F x {z) = zF (z). 

These two equations allow us to express Fq and Fx, and hence G\, in terms of F(z, 1): 

Gx(z) = (l-z)F(z,l)-l, 

This equation relates the two unknown functions of Eq. (|15|) . We replace G\{z) by the above 
expression in (|15j) . so that only one unknown function, namely F(z, 1), is left. The kernel 
method finally gives: 



, N 3 - 3z 2 - 2z 3 - 1 + z)y/l - 2z - 3z 2 , „ 2 o o ^ 4 in c ^. 

F(z, l = — ^ 4 jt = l + z + 2z 2 + 3z 3 + 6z 4 + I2z 5 + 0(« 

v ' 2(1 - z - z 2 + z 3 + z 4 ) 



Example 12. 

Let us choose A = {1}, B = {2} et C = {2}. The rewriting rule is now: 

(fc)-(0)(l)(3)(4)(5)...(fc-3)(fc-l)(fc + l). 
The functional equation reads 
u 2 (l - u) + zu 2 - zu 3 (l - u) + z(l - u) F(z, u) 



□ 
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= u 2 (l -u) + zu 2 F(z, 1) + z(l - u) \F Q {z) + uFi(z)} - zu 4 (l - u)G 2 (z), (16) 

with G2(z) = F(z,l) — Fq(z) — F\(z) — F 2 (z) — F±{z). Only three roots, uq,ui,u 2 can be 
substituted for u in the kernel, while the right-hand side of the equation contains four unknown 
functions, F(z, 1), F (z), F\{z) and G 2 (z). Writing (fTl)l for j = 0, 1 and 2 yields 

F (z) =l + z[F(z,l)-F (z)-F 2 (z)}, 

F 1 (z) = zF (z) + z [F(z, 1) - F (z) - F^z) - F 3 (z)} , 

F 2 (z) =zF 1 (z). 

The second equation is not of much use but, by combining the first and third one, we find 

l + z[F(z,l)-zF 1 (z)] 



Fo(z) 



l + z 



Replacing F (z) by this expression in ()16j) gives: 
u 2 {\ -u) + zu 2 - zu 3 (l - u) + z(l - u)l F(z, u) = u 2 (l-u) + 



z(l 



l + z 



+ zF(z,l) 



zl 



l + z 



+ z(l-u)F 1 {z) 



l + z 



zu\l-u)G 2 (z). (17) 



We are left with three unknown functions, related by three linear equations obtained by 
cancelling the kernel. Solving these equations would give F(z, 1) as an enormous rational 
function of z, uq,u\ and u 2 , symmetric in the Uj. This implies that F(z, 1) can also be 
written as a rational function of z and v = v\, the fourth and last root of the kernel. In 
particular, F(z, 1) is algebraic of degree at most 4. 

In order to obtain directly an expression of F(z, 1) in terms of z and v, we can proceed 
as follows. Let R'(z,u) denote the right-hand side of Eq. (|T7|). Then R'(z,u) is a polynomial 
in u of degree 5, and three of its roots are uo, ui,u 2 . Consequently, as a polynomial in u, the 
kernel K(z,u) divides (u — v)R'(z,u). 

Let us evaluate (u — v)R'(z,u) modulo K(z,u): we obtain a polynomial of degree 3 in 
n, whose coefficients depend on z, v, F(z, 1), F\ (z) and G 2 (z). This polynomial has to be 
zero: this gives a system of four (dependent) equations relating the three unknown functions 
F(z, l),F\(z) and G 2 (z). Solving the first three of these equations yields 

1 + z + z 2 - (z + l)zv + (z+ l)zv 2 



F{z,l) 



z 2 v 3 



1 — z 2 — z{l — z 2 )v + Z 3 V 3 
1 + z + 2z 2 + 3z 3 + 6z 4 + llz 5 + 23z 6 + 47z 7 + 101x 8 + 0(z 9 ) 



Eliminating v between this expression and K(z, v ) = gives a quartic equation satisfied by 
F(z,l). ' □ 



4 Transcendental systems 
4.1 Transcendence 

The radius of convergence of an algebraic series is always positive. Hence, one possible reason 
for a system to give a transcendental series is the fact that its coefficients grow too fast, so 
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that its radius of convergence is zero. This is the case for System (2.h), as proved by the 
following proposition. 

Proposition 6 Let b be a nonnegative integer. For k > 1, let m(k) = \{i : ei(k) > k — b}\. 
Assume that: 

1. for all k, there exists a forward jump from k (i.e., &i(k) > k for some i ), 

2. the sequence (ra(k))k is nondecreasing and tends to infinity. 

Then the (ordinary) generating function of the system has radius of convergence 0. 

Proof. Let so be the axiom of the system. Let us denote by h n the product m(so + b)m(so + 
2b) ■ ■ ■ m(sQ + nb). Let us prove that the generating tree contains at least h n nodes at level 
n(b + 1). At level nb, take a node v labeled k, with k > sq + nb. Such a node exists thanks to 
the first assumption. By definition of m(k), this node v has m(k) sons whose label is at least 
k — b. As m is non decreasing, v has at least m(so + nb) sons of label at least so + (n — 1)6. 
Iterating this procedure shows that, at level nb + i, at least m(so + (n — i + 1)6) • • • m(so + nb) 
descendants of v have a label larger than or equal to so + (n — i)b, for < i < n. In particular, 
for i = n, we obtain at level n(6 + 1) at least h n descendants of v whose label is at least sq. 

Hence f n (p+i) > h n . But as h n /h n -x = m(so + nb) goes to infinity with n, the series 
J2 n h n z n<yb+1 ^ has radius of convergence 0, and the same is true for F(z) = J2 n fnZ n - ■ 

In particular, this proposition implies that the generating function of any ECO-system in 
which the length of backward jumps is bounded has radius of convergence 0. Many examples 
of this type will be given in the next subsection, in which we shall study whether the corre- 
sponding generating function is holonomic or not. The following example, in which backward 
jumps are not bounded, was suggested by Nantel Bergeron. 



Example 13. A fake factorial walk 

Consider the system with axiom (1) and rewriting rules {(k) ~» (2)(4) • • • (2k)}. Proposition 
applies with 6 = and m(k) = 1 + L^/2j- Note that the radius of convergence of F(z) is 
zero although all the functions ej are bounded, and indeed constant: ei(k) = 2i for all k > i. 
The series F(z) is of course transcendental. Note, however, that F(z, u) satisfies a functional 
equation that is at first sight reminiscent of the equations studied in Section |3J 

r,/ x ^2 F(z,l)-F(z,u 2 ) 

F{Z, u) = U + ZU 7; • 

1 — IF 

□ 

The following example shows that Proposition |H1 is not far from optimal: an ECO-system 
in which all functions ej grow linearly can have a finite radius of convergence. 



Example 14. 

The system with axiom (1) and rules (k) (\k/2~\) k ~ 1 (k + 1) leads to a generating function 
with a positive radius of convergence. 

Let us start from the recursion defining the numbers f n ,k- We have _fo,i = 1 arid for n > 1, 

fn+l,k = /n,fc-l + (2fc — l)/ n ,2fc + (2fc — 2)f n2 k-l- 
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n k 
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3 
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3 


9 


7 


3 





1 


5 


1 





3 


7 


9 


3 



Table 1: The numbers f n ^ and g n ,k- Observe the convergence of the coefficients. 



The largest label occurring at level n in the tree is n + 1. Let us introduce the numbers 
9n,k = fn,n-k+i, for k < n. The above recursion can be rewritten as: 

9n+i,k = 9n,k + (2n - 2/c + 3)g ni2 fc-n-3 + (2n - 2fc + 2)g n ^k-n-2- (18) 

We have g nj fc = for k < 0. Hence Eq. (fTH)) implies that for k > 0, the sequence (g nj k)n is 
nondecreasing and reaches a constant value g(k) as soon as n > 2k — 1 (see Table 
Going back to the number f n of nodes at level n, we have 

n n 

fn = J2g n ,k < ^2g{k). 

k=0 k=0 

But 

n -.n 
n>0 fc=0 fe=0 

and hence it suffices to prove that the generating function for the numbers has a finite 
radius of convergence, that is, that these numbers grow at most exponentially. 
Writing (|18|) for n + 1 = 2k — i, for 1 < i < k, we obtain: 

92fc-i,fc = 92k-i-i,k + ( 2& - 2i + l)52fe-i-i,i-2 + (2fe - 2i)g 2 k-i-i,i-i ■ 
Iterating this formula for 7 between 1 and yields 

k 

g{k) = g2k-X,k = [( 2fc ~ 2i + 1 )fi'2fe-i-l,i-2 + (2fe - 22)52fc-i-l,i-l] 
i=l 

k k-2 

< J2 i( 2k ~ 2i + - 2 ) + ( 2k ~ 2i )a{i - i)] = E( 4A: " 4i " 5 )^(0- 

i=l i=0 

This inequality, together with the fact that g(0) = 1, implies that for all k > 0, ^(fc) < g(k), 
where the sequence g(k) is defined by g(0) = 1 and g{k) = J2i=o(^k — \i — 5)g(i) for k > 0. 
But the series J2kd( k ) zk ^ s rational, equal to (1 — z) 2 /(l — 2z — 2z 2 — z 3 ), and has a finite 
radius of convergence. Consequently, the numbers g(k) and g(k) grow at most exponentially. 
□ 

Algebraic generating functions are strongly constrained in their algebraic structure (by 
a polynomial equation) as well as in their analytic structure (in terms of singularities and 
asymptotic behaviour). In particular, they have a finite number of singularities, which are 
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algebraic numbers, and they admit local asymptotic expansions that involve only rational 
exponents. A contrario, a generating function that has infinitely many singularities (e.g., a 
natural boundary) or that involves a transcendental element (e.g., a logarithm) in a local 
asymptotic expansion is by necessity transcendental; see |16| for a discussion of such tran- 
scendence criteria. In the case of generating trees, this means that the presence of a condition 
involving a transcendental element is expected to lead to a transcendental generating function. 
This is the case in the following example. 

Example 15. A Fredholm system 

We examine System (2.e), in which the rules are irregular at powers of 2: 

(s ) = (2), (k)^(2) k ' 2 (3-[3p:k = 2P])(k + l), k>2. 

This example will involve the Fredholm series h(z) := J2 p >i z 2 \ which is well-known to admit 
the unit circle as a natural boundary. (This can be seen by way of the functional equation 
h(z) = z 2 + h(z 2 ), from which there results that h(z) is infinite at all iterated square-roots of 
unity.) According to Eq. ((21), we have, for k > 3, F k (z) = zFk-i(z), so that 

F k (z) = z k ~ 3 F 3 (z) for k > 3. 

Now, writing Eq. (J2J) for k = 2 gives 

F 2 (z) = l + zJ2(k-2)F k (z) + zJ2F2p(z) 

k>3 p>l 



l + zF 2 (z)+F 3 (z) (—L—+ h[ '- ) 



2 



For k = 3, we obtain: 



(l-zf z 



F 3 (z) = zF 2 (z) + z J2 F k( z ) 

fc>3, k^2f 

1 h(z) 



zF 2 (z)+F 3 (z) 



1 — z z 2 

Solving for F 2 (z) and F 3 (z), then summing (F(z) = F 2 (z) + F 3 (z)/(1 — z)), we obtain: 

F(z) = Z)2H ^} j = 1 + 2z + 5z 2 + Uz 3 + 39z 4 + 108z 5 + 0(x 6 ). 

(1 — 2z)(l — z) z h(z) — z* 

The functions h(z) and F(z) are rationally related, so that F(z) is itself transcendental. The 
series h has radius 1, but the denominator of F vanishes before z reaches 1 - actually, before 
z reaches 1/2. Hence the radius of F is the smallest root of its denominator. Its value is 
easily determined numerically and found to be about 0.360102. □ 



20 



4.2 Holonomy 



In the transcendental case, one can also discuss the holonomic character of the generating 
function F(z). 

A series is said to be holonomic, or D-finite |25| . if it satisfies a linear differential equation 
with polynomial coefficients in z. Equivalently, its coefficients f n satisfy a linear recurrence 
relation with polynomial coefficients in n. Consequently, given a sequence f n , the ordinary 
generating function J2 n f n z n is holonomic if and only if the exponential generating function 
J2 n fn,z n /n\ is holonomic. The set of holonomic series has nice closure properties: the sum or 
product of two of them is still holonomic, and the substitution of an algebraic series into an 
holonomic one gives an holonomic series. Holonomic series include algebraic series, and have 
a finite number of singularities. This implies that Example llf>| for which F{z) has a natural 
boundary, is not holonomic. 

We study below five ECO-systems that, at first sight, do not look to be very different. 
In particular, for each of them, forward and backward jumps are bounded. Consequently, 
Proposition El implies that the corresponding ordinary generating function has radius of con- 
vergence zero. However, we shall see that the first three systems have an holonomic generating 
function, while the last two have not. We have no general criterion that would allow us to 
distinguish between systems leading to holonomic generating functions and those leading to 
nonholonomic generating functions. 

Among the systems with bounded jumps, those for which ej(fc) — k belongs to {—1,0, 1} 
for all i < k have a nice property: the generating function for the corresponding excursions 
(walks starting and ending at level 0) can be written as the following continued fraction |15j : 



oxCqz 2 
l-b z 



-i j a 2 c 1 z 2 
1 — b\z 



a 3 c 2 z 2 
l-b 2 z 



where the coefficients a^jb^ and are the multiplicities appearing in the rules, which read 
(k)^> (k - l) ak (k) bk (k + l) Cfe . 

Example 16. Arrangements 

The system (k) {k){k + l) fc_1 with axiom (so) = (2) generates a sequence that starts with 
1, 2, 5, 16, 65, 326 (M1497). It is not hard to see that the triangular array f n< k+2 is given by 
the arrangement numbers k\(^j, so that the exponential generating function (EGF) of the 
sequence is 

F(Z,U)= J2 fn,kU k —r = ~ • 

^TT^o n - l — uz 

n>0,k>2 

This system satisfies the conditions of Proposition El with b = and m(k) = k. Accordingly, 
one has f n ~ enl, so that the ordinary generating function F(z) has radius of convergence 
and cannot be algebraic. However, F{z, 1) = e 2 /(l — z) is holonomic, and so is F(z). □ 



Example 17. Involutions and Hermite polynomials 

The system (k) ~> (k — l) k ^ 1 (k + 1) with axiom (so) = (1) generates a sequence that starts 
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with 1,1,2,4,10,26,76 (M1221). These numbers count involutions: more precisely, one 
easily derives from the recursion satisfied by the coefficients f n> i- that f n> k is the number of 
involutions on n points, k — 1 of which are fixed. Proposition |H1 applies with 6 = 1 and 
m(k) = k. 

The corresponding EGF is 



z 2 



F(z,u) = Y fn,ku k — =uexp i zu + — , (19) 

n>0,fc>l U - V 1 ) 

and its value at u = 1 is holonomic. 

The polynomials f n (u) = J2k fn,kU k counting involutions on n points are in fact closely 
related to the Hermite polynomials, defined by: 

Hnix)*— = exp [st-M. 

n>0 H - \ Z / 

Indeed, comparing the above identity with (|19|) shows that f n (u) = u i n H n (—iu). □ 



Example 18. Partial permutations and Laguerre polynomials 

The rewriting rule (k) (k + \) k ~ l {k + 2), taken with axiom (2), generates a sequence that 
starts with 1,2,7,34,209,... (M1795). From the recursion satisfied by the coefficients f n k, 
we derive that f n ,n+k is the number of partial injections of {1,2, ... ,n} into itself in which 
k — 2 points are unmatched. From this, we obtain: 

F(z, «) = exp (-^-) =u^ Ln(-u)^f- 

1 — uz \1 — uz ^— ' n! 

\ / n>0 

where L n {u) is the nth Laguerre polynomial. Again, F(z, 1) is holonomic. □ 
The next two systems, as announced, lead to nonholonomic generating functions. 



Example 19. Set partitions and Stirling polynomials 

Let us consider the system [(l),(fc) ~> (k) k ~ 1 (k + 1)]. From the recursion satisfied by the 
coefficients f Ui k, we derive that f n ,k+i is equal to the Stirling number of the second kind {^}, 
which counts partitions of n objects into k nonempty subsets. The corresponding EGF is 

F(z, u) = u exp (ii(exp z — 1)) . 

At u = 1, this generating function specializes to 

F(z, 1) = exp(exp(z) - 1)) = £ B n ^ = 1 + z + 2^ + 5^ + 15^ + 52^ + 203^ + . . . 

n>0 

This is the exponential generating function of the Bell numbers (M1484). It is known that 
logB n = n log n — n log log n + 0(n) (see [20 ), and this cannot be the asymptotic behaviour of 
the logarithm of the coefficients of an holonomic series (see |29) for admissible types). Hence, 
F(z, 1), as well as F(z, 1), is nonholonomic. □ 
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Example 20. Bessel numbers 

We study the system with axiom (2) and rewriting rules 



(2)~>(2)(3), (k)^ (k-l)(k) k - 2 (k + l), k>3. (20) 

We shift the labels by 2 to obtain a walk model with axiom (0) and rules 

(0)~>(0)(1), (k) ^ (k — l)(k) k (k + 1), k>l. 

The corresponding bivariate generating function F(z, u) satisfies the functional differential 
equation 

F(z, u) ( 1 - z(u + n' 1 )) =l + z(l- u' 1 )F(z, 0) + zu—(z, it), 

V / OU 

which is certainly not obvious to solve. However, as observed in ^5], it is easy to obtain a 
continued fraction expansion of the excursion generating function: 

F(z, 0) = l+z+2z 2 +iz 3 +9z 4 +- ■ ■ = 1 ' 



z 



2 ~\-z-z 2 B{ z y 



z 2 



l-2z — 



1 - 3z - '■■ 

where B(z) = E n B n zn = 1 + z + 2z 2 + 5z 3 + Uz 4 + 43z 5 + U3z 6 H is the generating 

function of Bessel numbers (M1462) and counts non-overlapping partitions As F(z,0) 
itself, the series B{z) has radius of convergence zero. The fast increase of -B* entails 

[z n ]F(z,0) ~ Bl_ 2 . 

From ^7j, we know that log-B* = nlogn — n log log n + 0(n). Again, this prevents F(z,0) 
from being holonomic. 

In order to prove that F(z, 1) itself is nonholonomic, we are going to prove that its 
coefficients f n have the same asymptotic behaviour as the coefficients of F(z,0). Clearly, 

[z n ]F(z,0) =/„,0<£/n,fc = /n. 

k 

To find an upper bound for /„, we compare the system (|20|) (denoted Si below) to the system 
£2 with axiom (2) and rule (k) ~> (k) k ^ 1 (k + 1). This system generates a tree with counting 
sequence g n . The form of the rules implies that the (unlabeled) tree associated with Si is a 
subtree of the tree associated with £2. Hence f n < g n . Comparing S2 to the system studied 
in the previous example shows that g n is the Bell number B n+ i, the logarithm of which is also 
known to be nlogn — n log log n + 0(n) (see |20|). Hence log/ n = nlogn — n log log n + 0(n), 
and this prevents the series F(z, 1) from being holonomic. □ 
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Axiom 



System 



Name 



Id. 



Generating Function 



(1) 

(2) 
(3) 



Rational OGF 

(ft) ^ (fc) fc - 1 ((/c mod 2) + 1) 
(ft) (2) k - 1 (k + l) 
(k) ~> (2) k ~ 1 (k + 1) 



Ex. Fibonacci 
Ex. 0] even Fibonacci 
Ex. 0] odd Fibonacci 



M0692 
M1439 
M2741 



OGF 

i 

1-z-z 2 

1-z 
l-3z+z 2 
1 

l-3z+z 2 



(1) 

(2) 
(3) 
(4) 
(to) 

(3) 
(4) 
(to) 



Algebraic OGF 



(k) 

(k) 

(k) 

(k) 

(k) 

(ft) 
(k) 
(k) 



(1) 
(2) 
(3) 
(4) 
(to) 

(3) 
(4) 
(to) 



(fe-l)(fc+l) 



(k)(k + 
(k)(k + 
(k)(k + 

•(A)(fcH 

(ft + 2) 
(ft + 3) 
■ (ft + m 



1) 
I) 2 

1)3 
l) r ' 



1) 



Ex. [7| Motzkin numbers 
Ex. H2 Catalan numbers 
Ex. Schroder numbers 



Ex. OH Ternary trees 

Ex. El Dissections ol a polygon 

Ex. TO-ary trees 



M1184 
M1459 
M2898 
M3556 



M2926 
M3587 



OGF 

l-z-\/l-2z-3z 2 



l-22-Vl-4z 
2P 



l-3z-Vl~6z+z 2 
IP 



l-4z-Vl-8z+4z 2 
6z 2 



1 — mz — -^/l — 2mz+(m — 2) 2 z 2 
2(m-l)z 2 

F = (1 + zF) 3 
F= (1 + zF) 4 
F = (1 + zF) m 



(1) 

(2) 

(1) 
(2) 
(2) 
(2) 



Holonomic 
transcendental OGF 

(ft) (k + l) k 
(k)^ (k)(k + l) k - 1 



ik-l 



(ft. 
- l {k- 



(fe) 
(k) 

(*)■ 

(fc)~» (ft- l) fe - 2 (fc + l) 



(ft. 
(fc- 

(fc-l) fe - 2 (fc)(fc + l) 

fc-2/ 



Permutations 

Ex. El Arrangements 

Ex. El Involutions 

Ex. El Partial permutations 

Switchboard problem 

Bicolorcd involutions 



M1675 
M1497 
M1221 
M1795 
M1461 
M1648 



EGF 

1/(1 -z) 
e z /(l - z) 

z+Iz 2 

/(l- 

3 2z+iz 2 
„2z+z 2 



e 

,z/(l-z) 



(1) 

(2) 
(2) 



Nonholonomic OGF 

(k) k - 1 (k + 1) 
(fc)-> (fc) fc - 2 (fc + l) 2 
(fc)-> (A - l)(fc) fe - 2 (fc + l) 



Ex. El Bell numbers 
Bicolored partitions 
Ex. E| Bessel numbers 



M1484 
M1662 
M1462 



EGF 

=2(e a -l) 



Table 2: Some ECO-systems of combinatorial interest. 



A small catalog of ECO-systems 

To conclude, we present in Table 2 a small catalog of ECO-systems that lead to sequences 
of combinatorial interest. Several examples are detailed in the paper; others are due to West 
|27[ I28j or Barcucci, Del Lungo, Pergola, Pinzani [U El El E| > or are folklore. Each of them is 
an instance of application of our criteria. 

Acknowledgements. We thank Elisa Pergola and Renzo Pinzani who presented us the 
problem we deal with in this paper. We are also very grateful for helpful discussions with 
Jean-Paul Allouche. 
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